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INTRODUCTION 


Because  magnetoplasmadynamic  (MPD)  thrusters  are  inherently  low  thrust  devices, 
operation  for  hundreds  to  thousands  of  hours  will  be  required  to  impart  useful  levels  of  total 
impulse.  Experience  at  Princeton1-2  and  in  other  laboratories^4  indicates  that  the  cathode  of  these 
devices  suffers  the  most  severe  damage  in  the  hostile  environment  of  the  high  current  discharge, 
therefore  representing  the  life-limiting  component.  To  explore  the  physical  mechanisms 
responsible  for  cathode  degradation,  a  new  diagnostic  tool  has  been  developed,  the  surface  layer 
activation  (SLA)  technique,  to  monitor  erosive  loss.  This  method  was  chosen  in  the  Phase  I 
portion  of  this  contract  from  several  alternatives  because  it  provides  highly  accurate,  in-situ 
measurements  of  erosive  mass  loss  and  can  be  readily  applied  to  a  wide  variety  of  materials  and 
environments.  Concurrent  programs  of  erosion  measurement  on  experimental  MPD  thrusters 
using  this  technique  and  supporting  analytical  modelling  of  the  erosion  process  are  now  being 
pursued. 

The  SLA  diagnostic  technique  reues  upon  the  production  of  a  radioactive  tracer  in  a  thin  layer 
near  the  cathode  surface  by  nuclear  activation.  Monitoring  changes  in  the  activity  level  during 
thruster  operation  provides  a  direct  measure  of  the  amount  of  the  activated  layer  removed  by  the 
destructive  plasma-surface  interaction.  Partial  results  of  preliminary  tests  of  the  technique  have 
already  been  reported1 .  The  purpose  of  this  paper  is  to  provide  the  completed  set  of  preliminary 
data,  discuss  substantial  improvements  w  hich  make  surface  layer  activation  a  practical  tool,  and 
describe  the  current  direction  of  the  experimental  program. 

In  the  interests  of  clarity  and  conciseness,  many  details  of  the  work  have  been  omitted.  A 
comprehensive  archival  literature  base  exists  for  this  work  in  the  form  of  monthly  technical 
progress  repons  which  are  available  from  the  Princeton  University  Engineering  Library.  The 
documents  for  this  reporting  period  have  been  bound  into  yearly  compendia  which  can  be  obtained 
by  requesting  MAE  Report  #1678  for  1985  and  MAE  Report  #1679  for  1986.  Alternatively, 
individual  monthly  reports  relating  to  specific  topics  car  be  obtained  directly  from  A.J.  Kelly. 

PRELIMINARY  TESTS  ON  A  HALF-SCALE  FLARED  ANODE  THRUSTER 
APPARATUS  AND  PROCEDURE 

The  coaxial  thruster  configuration  shown  in  Figure  1  was  chosen  for  the  second  test  of  the 
surface  layer  activation  technique  because  preliminary  performance  data5  indicated  that  this  design 
is  superior  to  the  benchmark  configuration  studied  in  the  first  test.  In  this  design,  a  high  speed 
solenoid  valve  controls  the  injection  of  argon  through  the  boron  nitride  backplate  via  an  annulus 
surrounding  the  central  2%  thoriated  tungsten  cathode  and  by  12  small  holes  arrayed  on  a  circle 
near  the  copper  anode.  The  thruster  was  mounted  in  a  fiberglass  vacuum  tank  1.8  m  in  diameter 
and  4.8  m  long,  which  was  typically  maintained  at  a  pressure  of  about  10'5  torn  One  msec 
duration  rectangular  current  pulses  were  supplied  by  a  3000  pF,  175  kJ  pulse-forming  network. 

The  activation  process  will  be  discussed  in  detail  in  the  next  section;  however,  a  brief 
description  of  the  activation  and  data  reduction  used  in  this  experiment  is  included  here  to 
demonstrate  several  deficiencies  which  have  subsequently  been  corrected.  Gamma  emitting 
radionuclide  tracers  were  produced  by  activating  small  spots  5  mm  in  diameter  on  the  tungsten 
cathode,  copper  anode,  and  boron  nitride  backplate  at  the  Princeton  Cyclotron  using  the  parameters 
listed  in  Table  1 .  These  spots  were  placed  at  positions  of  maximum  current  density5  where  the 
most  severe  erosion  presumably  occurs  The  resulting  activity  was  monitored  by  observing 
prominent  photopeaks  in  the  gamma  ray  spectrum,  which  was  collected  with  standard  radiation 
detection  equipment  mounted  outside  the  thruster.  Gamma  rays  absorbed  by  a  Bicron  5x5  inch 
Naim )  scintillator  crystal  produce  a  number  of  photons  proportional  to  the  gamma  energy.  The 
photons  are  detected  with  a  photomultiplier  tube,  which  produces  voltage  pulses  that  are  then  read 


and  displayed  as  the  number  of  detected  decay  events  (counts)  in  a  particular  energy  range 
(channel)  by  a  Canberra  Series  35  multichannel  analyzer.  The  composite  gamma  ray  spectrum  ot 
the  lab  background  and  the  three  isotopes  produced  in  the  activation  is  shown  in  F  igure  2 

The  activity  of  each  source  is  proportional  to  the  number  of  counts  in  the  photopeak  .which 
was  determined  by  subtracting  the  background  which  is  assumed  to  be  linear  and  summing  the 
remaining  counts,  as  shown  in  Figure  3.  To  relate  the  fraction  of  the  original  activity  remaining  to 
the  amount  of  material  removed  by  erosion,  the  distribution  of  the  radioactive  tracer  in  the  activated 
layer  must  be  known.  This  was  determined6  by  repetitively  lapping  thin  layers  from  an  irradiated 
sample  of  the  material  of  interest  and  determining  die  remaining  activity  to  produce  a  calibration 
curve  such  as  that  plotted  in  Figure  4. 

The  test  sequence  for  this  experiment  consisted  of  10,000  1  msec  long  discharges  at  an  argon 
mass  flow  rate  of  2.2  g/s  and  a  current  level  of  9.6  kA  (below  the  onset  current),  followed  by 
4000,  1  msec  discharges  at  0.75  g/s  and  12  kA  (above  onset).  The  activities  of  the  three  isotopes 
were  measured  periodically  during  this  sequence  to  determine  the  cumulative  material  losses. 

RESULTS 

No  anode  or  backplate  erosion  was  detected  for  either  operating  condition.  The  cumulative 
cathode  mass  loss  per  unit  surface  area  is  plotted  in  Figure  5  as  a  function  of  the  number  of 
discharge*  frir  both  operating  conditions.  Figure  6  displays  the  same  data  plotted  against  the  total 
charge  transferred  through  the  cathode;  however,  because  the  eroded  mass  is  more  likely  a  function 
of  the  local  current  density  in  the  activated  spot,  which  may  have  differed  for  the  two  operating 
conditions,  comparisons  of  the  two  sets  may  not  be  valid.  The  most  striking  characteristic  of  the 
data  is  that  regions  of  strongly  linear  behavior  are  separated  by  regions  of  discontinuous  mass  loss, 
a  behavior  not  found  in  the  first  application  of  the  technique1^.  As  shown  in  the  figures,  the 
continuous  regions  are  well  represented  by  straight  line  fits  which  yield  the  erosion  rates  listed  in 
Table  2,  which  are  similar  to  weight  loss  measurements  made  at  Stuttgart7  and  Tokyo3. 

Although  no  surface  distortion  in  the  activated  spot  was  visible  on  a  macroscopic  or 
microscopic  scale  prior  to  the  test,  after  the  14,000  discharges  the  spot  was  cleaily  visible  as  a 
shallow  depression  with  a  flat  floor  and  a  nm  of  tungsten  that  apparently  melted  and  peeled  a\>  ay 
from  the  center.  A  similar  blister  was  found  on  the  cathode  used  in  the  first  test  of  the  technique. 
Inspection  of  the  cathode  surface  under  an  optical  microscope  revealed  a  number  of  pits  or  craters 
which  were  densely  parked  and  very  distinct  at  the  tip,  less  dense  ard  distinguishable  in  the  center 
and  very  indistinct  at  the  base,  where  large  scale  surface  melting  and  flowing  had  occurred.  Under 
moderate  magnification,  small  metallic  spheres,  apparently  tungsten,  were  visible  on  the  backplate 
near  the  annular  injection  port  which  surrounds  the  cathode. 

Much  smaller  scale  cratering  was  evident  under  a  scanning  electron  microscope,  typified  by 
the  photomicrograph  in  Figure  7.  The  typical  crater  diameter  was  1  to  5  microns  at  the  tip  and  the 
base,  and  1-2  microns  in  the  center.  No  qualitative  differences  or  general  differences  in  scale  could 
be  discerned  on  a  microscopic  level  between  the  irradiated  spot  and  regions  surrounding  it  O- 
the  larger  pits  visible  in  the  optical  microscope  was  examined  with  the  SFM,  and  appeared 
as  a  depression  60  70  microns  in  diameter  thickiy  overlaid  with  the  snider  craters,  ou’  m  .ion 
wide  cracks  in  the  surface  were  visible  in  ali  regions. 

DISCUSSION 

The  appearance  of  the  activated  spot  after  the  experiment  strongly  suggests  that  atypical 
surface  damage  occurred,  probably  due  to  material  property  changes  induced  by  the  high  energy- 
ion  bombardment  which  produced  the  radioactive  tracer.  The  fraction  of  tungsten  transmuted  to 
,850s  was  0(10  17).  so  its  presence  as  an  impurity  could  not  be  responsible.  A  more  likely 
explanation  is  tha*  radiation  damage  from  tiie  23  MeV  a-particles  used  tor  the  activation  caused  the 
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rupture  of  the  surface  during  the  experiment  Gross  loss  of  material  from  the  peeled  back  edges  of 
the  region  may  have  been  responsible  foi  the  discontinuous  jumps  in  the  erosion  data.  Concern 
that  the  measurement  process  ts  significantly  disturbing  the  measured  phenomenon  motivated  a 
critical  examination  of  the  activation  process,  which  culminated  in  the  defi.ition  of  criteria 
governing  the  choice  of  an  activation  scheme  capable  of  providing  honest  estimates  of  surface 
degradation.  The  results  of  that  study  are  summarized  in  the  next  section. 

The  validity  of  the  erosion  rates  obtained  from  the  continuous  portions  of  the  data  are 
defensible  on  the  basis  of  the  qualitative  similarity  between  the  microscopic  erosion  structures 
found  in  the  activated  spot  and  the  undisturbed  regions  around  it  and  the  reasonable  agreement  with 
other  measurements.  Assuming  similar  ero-ion  rates  over  the  entire  cathode  yields  an  erosion  rate 
of  about  1  fig/C,  or  approximately  0.01  g/s  a:  10  kA,  which  is  clearly  unacceptably  high  for  a 
practical  device. 

On  cold  cathode  surtaces  current  continuity  is  maintained  primarily  through  a  number  of 
small  emission  sites  where  the  local  temperature  is  extremely  high,  liberating  electrons  from  the 
surface  through  a  combination  of  thermal  and  field  emission.  The  photomicrographs  of  the 
cathode  surface  confirm  that  the  erosion  processes  occur  primarily  in  localized  microspots 
corresponding  to  the  emission  sites  where  temperatures  exceed  the  melting  temperature.  The  photo¬ 
micrographs  showing  molten  tungsten  splashed  from  these  sites  and  the  tungsten  droplets  found 
on  the  backptiio.  ndicam  that  droplet  erosion  may  be  a  significant  component  of  the  mass  flux  from 
the  cathode  surl^o,  oiu,  vgh  the  high  local  temperatures  also  imply  high  evaporation  rates.  The 
appearance  of  larger  pits  aroun  *  ’00  microns  in  diameter  and  the  smaller  craters  1-5  microns  in 
diameter  suggests  that  two  types  of  arc  attachment  may  occur.  This  is  consistent  with 
Rakhovskii's  observations  of  a  vacuum  '  scharge  using  high  speed  photographic  techniques8.  He 
found  two  types  of  iuminous  spots-fast  moving  spots  that  left  small  craters  and  slower  ones  that 
tended  to  cluster  and  cause  more  extensive  damage. 

Evidence  exists  thui  when  the  bulk  cathode  surface  temperature  becomes  sufficiently  high, 
diffuse  thermionic  emission  provides  sufficient  current  and  the  local  melting  associated  with  the 
micros/Ot  mode  of  emission  does  not  occur,  allowing  an  erosion  rate  several  orders  of  magnitude 
lower  than  that  experienced  with  cold  cathode  bulk  temperatures7.  This  provides  some  hope  that 
steady  state  thrusters  in  which  the  cathode  is  heated  by  the  discharge  may  have  acceptable  cathode 
lifetimes.  The  next  stage  of  experimentation  described  in  the  final  section  of  this  paper  is  designed 
to  explore  this  regime  of  cathode  operation. 

ANALYSIS  OF  THE  SURFACE  LAYER  ACTIVATION  TECHNIQUE 

Two  critical  operations  compose  the  the  SLA  mass  loss  measurement  technique.  First,  a 
suitable  gamma-emitting  tracer  must  be  produced  in  a  thin  surface  layer,  second,  the  change  in 
aenvity  of  this  source  as  it  erodes  must  be  determined.  The  activation  and  spectrum  analysis 
fundamentals  will  be  introduced  next  to  establish  the  processes  and  terminology  used  subsequently 
in  defining  criteria  for  successful  application  of  the  rechnique.  Finally,  a  new  activation  scheme 
will  be  presented  which  more  satisfactorily  meets  these  criteria  than  that  used  in  the  preliminary 
tests,  and  which  is  conceptually  capable  of  providing  a  direct  measure  of  the  extent  of  surface 
pitting. 

PRODUCTION  OF  THE  RADIOACTIVE  TRACER 


The  key  to  the  SLA  technique  is  the  depth  calibration  curve,  which  relates  the  observable 
aenvity  decline  to  the  desired  measurement  of  mass  loss.  As  intimated  earlier,  this  relationship  is 
dependent  on  the  distribution  of  the  radionuclide  beneath  the  surface.  Because  this  technique  relies 
on  a  nuclear  reaction  in  the  target  surface  excited  by  a  high  energy  ion  beam,  the  resulting 
distribution  is  dependent  on  the  physics  of  that  interaction.  The  density  of  the  radioactive  atoms  na 
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is  the  number  of  target  atoms  that  undergo  the  particular  nuclear  react;' >'t  during  bomba/dor  -• 

na(x.v.z)  J  ^  n,,<x.y fT((F )  d>fx.y,7.t:  di 

where  nt(  is  the  density  of  isotope  i  in  the  target  material,  rr,  is  the  cross  section  fur  reach  or.*, 
ixrtween  the  beam  particle  of  energy  E(x,y,z)  and  target  isotope  1.  <?  is  the  flux  ot  beam  pane 
and  t  is  the  total  irradiation  time.  The  sum  is  taken  over  ail  reactions  between  die  bean,  panicle  a 
isotopes  in  the  target  which  produce  the  radioisotope.  The  x  ,  y  .  and  z-axes  form  an  orthogon.: 
coordinate  system  originating  at  the  intersection  of  the  beam  and  the  target  surface,  with  the  z-av 
directed  into  the  target  along  the  beam  axis.  Assuming  isotropic  target  density,  a  constant  ion 
energy  across  die  beam,  and  an  isotropic  and  time  independent  beam  density  (no  variation  m  dr 
incident  he  am  and  negligible  loss  of  particles  by  interactions  along  the  path  of  interest  in  'he  targ 
yields 


na(z)  --  £  ntl  o(Fd  <6  t 

Ot  represents  the  total  dose,  which  can  be  measured  during  activation.  Reaction 
cross-sections  are  often  available  from  measurements  or  theoretical  calculi  tr ‘ns.  lor  ;;v.ur'  c. 
Figure  8  shows  the  cross-section  for  the  ^Cuip.t; d(’^Zr.  reaction'  !  Since  the  cross  see  ten 
junction  of  the  beam  particle  energy,  ifie  z -dependence  of  the  energy  in  the  targe'  must  lx* 
determined  to  calculate  density  profiles  along  z. 

At  enei gies  above  about  1  MeV,  the  incident  panicle  'uses  energy  primarily  through 
interaction  with  the  electrons  of  die  target  material.  This  energy  loss  can  be  reliably  '.•aScu'nvd  : 
within  about  )%  ,2‘14.  The  energy  loss  can  be  integrated  along  the  actual  particle  path  to  obtain 
energy,  but  the  projection  of  the  actual  path  onto  the  z- axis  must  be  determined  for  on;  pi.rp-  . 
Fortunately,  large  angle  scattering  is  rare  ir.  encounters  with  electron-;  at  high  energies,  :o  a rim 
MeV'  the  path  of  the  incident  ion  is  very  nearly  straight’^  Below  1  MeV.  correc  tion  .  for 
interactions,  with  the  target  nuclei  and  for  deviations  from  a  straight  path  must  be  made.  Since  d 
nuclear  reactions  of  interest  have  threshold  energi  -s  above  1  Me\  however,  this  is  urinec-wca" 
Figure  tJ  dtsplavs  the  calculated  proton  <*  :v*rgv  in  copper  assuming  the  *  M  Me'  Incident  e:>e*. 
used  in  the  copjier  activali-  1  lie  extent  >r  the  activated  layer  is  deter:. •-med  by  the  point  in  the 
target  where  the  incident  particle  energy  drops  to  the  reaction  threshold  energy,  as  illustrated  in 
Figure  9. 

The  der":*y  referenced  •’.*  the  rj-a  .  roonal  to  :he  surface  and  directed  inv.  :;rd  i*-  ;mph 

related  *r>  !r f *'{.■?  }  \  !*»•«  *  '■*•  *rj‘T  •- •*  jfv»  nr ■  "  0  |v!\ r**n  (ht-  j* :  **•, 
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activated  area,  which  is  in  fact  the  process  used  to  generate  empirical  depth  calibration  curves. 

This  calibration  relation  is  expressed  theoretically  as 

r,  . 

^(PjTi'A)  =  1  j  nairj)  dq  )  /( {  na(T|)  drj  )  (IS) 

o  o 

where  £(ptr)'A)  is  the  fractional  activity  remaining  after  a  layer  of  thickness  rj'  and  area  A  has  been 
removed  and  q0  is  the  depth  of  the  activated  layer  measured  perpendicular  to  the  surface.  Of 
course,  for  data  analysis  the  inverse  relation  is  required.  Figure  10  demonstrates  the  agreement 
between  the  measured  depth  calibration  corves  and  the  curve  generated  by  integrating  the  calculated 
density  for  the  copper  activation.  The  dose  correspondence  between  the  modelling  and  experiment 
indicates  that  the  ultimate  utility  of  theoretically  generating  depth  calibration  curves  is  limited  only 
by  the  availability  of  appropriate  cross-section  data. 

ANALYSIS  OF  THE  RESULTING  GAMMA  SPECTRUM 

Two  methods  of  collecting  and  analyzing  activity  change  data  have  been  developed  The  first 
relies  on  a  comparison  of  the  gamma  spectrum  after  operation  of  the  thruster  with  the  initial 
spectrum,  which  yields  the  total  activity  change,  occurring  during  operation.  This  method  is  simple 
and  highly  accurate  if  the  correct  approach  is  used,  as  described  next.  The  second  method  involves 
monitoring  the  countrate  of  a  portion  of  the  spectral  region  during  operation,  and  while  it  is  not  as 
accurate,  it  allows  real-time,  time-resolved  mass  loss  measurement. 

The  spectrum  analysis  method  used  in  the  preliminary  tests  relied  on  a  linear  approximation 
for  the  background  spectrum,  an  assumption  that  is  in  general  not  justified,  particularly  when 
photopeaks  overlap,  as  the  7Be  and  ^5Os  peaks  do  in  Figure  2.  Not  the  least  worrisome  aspect 
of  this  approximation  is  that  the  error  cannot  be  estimated  reliably.  The  following  method  is  much 
more  sound  theoretically  and  has  been  tested  and  found  to  eliminate  virtually  all  of  the  uncertainty 
in  the  activity  measurement. 

The  assumptions  underlying  this  method  of  spectrum  resolution  are15; 

a)  The  sample  spectrum  is  produced  by  a  combination  of  known  isotopes. 

b)  The  response  function  of  the  detector  for  each  isotope  is  known  and  is  independent  of 
the  activity  level. 

c)  The  sample  spectrum  is  a  linear  combination  of  the  response  functions  of  the 
component  isotopes. 

d)  Each  component  has  a  different  spectrum,  all  of  which  are  linearly  independent. 

The  first  of  these  is  easily  met  in  this  experiment,  since  the  isotopes  produced  in  the  activation  are 
known.  The  response  function  required  in  (b)  can  be  determined  by  measuring  the  shape  of  the 
gamma  spectrum  of  an  isolated  source  with  the  same  detector- source  geometry  to  be  used  in  the 
experiment,  provided  the  count  rate  is  low  enough  to  avoid  coincidence  summing  errors,  (c)  is 
also  valid  i*-  the  count  rate  is  low  enough,  and  can  be  expressed  as 

Ni  =  £  xk  Aki  +  fi  (6) 

k 

where  Nj  is  the  number  of  counts  in  channel  i,  is  the  intensity  of  radioisotope  k  relative  to  the 
measured  response  function  (the  fractional  activity,  if  the  response  function  is  measured  from  the 
activated  component  prior  to  erosion),  Ay  is  the  number  of  counts  in  channel  i  in  the  response 
function  for  isotope  k,  and  £j  is  the  measurement  error.  The  background  can  be  considered  one  of 
the  source  spectra  or  measured  separately  and  subtracted  out  to  produce  pure  source  spectra. 
Assumption  (d)  must  be  met  by  a  careful  selection  of  isotopes. 
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A  linear  least  squares  program  is  •;  .eu  to  obtain  the  best  estimate  of  the  mten  ,ue>  .» j, 

1  he  emu  is  inversely  proportional  to  the  number  of  count*  >rt  each  ..  'nnao.  and  i  a-  -K  red.;, 
an  arbitrarily  low  value  bv  ipcrcasipp  the  counting  time  or  avetaginj.  over  a  numhe:  of  . 
i  'rrors  less  than  0..V*  have  been  obtained  routinely  with  ntitnvune  level  sown  rs!w 

Instead  of  analyzing  the  photopeaks  t<>  the  speemtm  c olleeted  after  opera  mu.  c  th-  tl  r. 
which  yields  the  time  integrated  mass  loss,  the  total  counts  in  a  particular  energy  range  can  b 
monitored  as  a  function  of  time  The  decrease  in  the  counting  rate  with  tune  is  then  a  direct 
measure  of  the  materia!  lost  during  that  time.  The  multichannel  analyzer  used  in  the  conventiot 
analysis  method  can  be  programmed  to  operate  in  this  mode  by  initially  setting  upper  and  low- 
level  discriminators  to  filter  out  ali  pulses  outside  a  specified  range  and  then  displaying  u: 
successive  channels  the  puises  counted  during  a  certain  dwell  time  The  pulse  height  range  err 
chosen  to  isolate  the  spectra!  region  with  the  most  intense  photopeaks,  which  will  maximize  "  • 
signal-to- noise  ratio,  and  the  dwell  time  per  channel  can  be  varied  from  microsecond'  to  hom 
give  an  appropriate  number  of  counts  per  channel  and  sufficient  time  resolution  for  the  crush  >: 
process  being  studied  Preliminary  estimates  confirm  that  modest  levels  of  activity  should  p 
data  with  a  fairly  low  level  of  uncertainty  for  time  resolution  down  to  aUuit  P.  1  >e<  >nj 

CRITKRI A  K)R  Ob'WT  Of-  \\  At  T  M!«'n  v  m  vij. 

1  he  following  pnrcuu  • :  ■  •< wsyieratn  re.e’i  f-r  .  ioi.-e  n*  nmvj'ziH  v  u 

characterized  by  the  beam  eucrgv  and  top.  •.  'he  ■}.  n;«.  -nd  macnttmli  of  ' t  •  •  c-o-  •  -  .!-• 

threshold  energv  for  th'-  reaction,  'he  pr  •  '  'h-  -*- u :  -.n  pro!  .  "...  ,:nd  f  <•  oo'lod  1 

spectrum  resolution,  as  explained  above 

L-High  v 


Ihe  yield  Y  can  be*  expressed  as  the  n umber  ot  radioisotop**  atoms  produced  pc  uni' 
deposited  beam  m.itcrud  m  !  is  g nvp  •. ■<  bv  the  "perpv  of  the  beam  .md  tbe  rv'<p.»vt,'de  o’  i 
cross-section-  i n ■  i •_  • 

>  |  2.  nt,  op, bid/ 

The  same  value  can  lx:  c  ptesst-d  in  to-  u  ,■  .n  t ■ . sal  teirns  as  the  number  of  niicTocuric-'  •>*  a 
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displace^  in;  T  i  i  •  a,  inru;  .  •  u;  p.  article  with  ;•  pii-.tv’ 

knock -on  target  ;»»  mi  -phii  ’he  *.;>>•  .  •  .*  fhiuc  p-vp'  def-  t  •<-•.  on  bint  durum  the 

diffusion  that  recurs  even  at  nxm  temper.’ v.i  .  ,  how  ever,  some  wiM  dust.:,  together  to  form 
extended  defect-  s-  •;**  :n  did'  cai'ni:  he  ;■  '  'up  ar  !  cavuies  vtl.r.li  arc  mu.  h.  mote  v.mic 

Tucrostn'ctiini!  change-  r-  cd  bv  the  f««  -uunc-n  c  extended  dctlv ..  .rod  u.c  segregation  oi 
impurities  c>r ->l!ov' i  pf-t ;  iis  w h i  h  cs  »u  ompanv  ?he  4  ondfuxit:  >"  ■  i>  cef c. .  •  v;v,  r-" 
in  macroscopic  matcnal  p-onert;  - ‘up-p*'  -  i  ’  !■  am  r-urrred  to  a  rshp.-,  r  tut-  11  i<:inv 

.'  pur  . cue:.-,  initially  reared  *--y  an  >,  v  r  can  bu  ■  donated  ••..  ;■■.  .:■■  dc  .  rib  .’  in  n,- 
levs  paragraphs  Hei  a.;  tt;c  h-tmarion  •’  ■  ven  !cj  uc'ccis  dcpctiu--  ten  the  •Jeu.iied  kittens--.  c1 

vacancy  and  interstitial  migration  and  interaction  with  ranks,  theoretical  prediction  of  the  »c*ulii 
micTostn-'  tup’!  and  micnx.  op  spositiomd  is  mu.'!1  more  difficult.  >o  thru  extensio,.  o'  th 

analysis  will  not  be  pursued.  It  should  he  noted  qualitatively,  however,  that  die  presence  of  th 
highly  insoluble  gr.  -  helium  n  the  damaged  lattice  at  u  ncentrations  as  low  as  tens  o'  ppm  can 
significantly  aggravate  cavity  tomation  anti  the  resulting  effects**’. 

The  ph”  bn  ’>u  si'  •  >>!  g’inu:  v  k-s,  c  !*  calculated  u-  ■;  '(■  fo; 


npka  -  ^  1  "t  ad 


(X) 

which  is  similar  to  that  used  to  calculate  the  density  oi'  activated  atoms.  Here  nt  is  the  density  of 
tui  tter  atoms  and  oj  is  the  cross-section  for  displacing  lattice  atoms.  The  cross-section  for  energy 
transfer  greater  than  that  required  for  displacement  can  he  approximated  for  energies  above  a  few 
keV  using  the  Rutherford  formula,  which  describes  Coulomb  interactions  between  the  bombarding 
and  target  nuclei’*1.  This  equation  is  inversely  proportional  to  incident  particle  energy,  reflecting 
the  fact  that  at  high  energies  most  of  the  energy  loss  is  due  to  interaction  with  the  target  atom 
electrons,  not  the  nuclei.  This  guarantees  that  most  of  the  radiation  damage  will  occur  near  the 
mean  depth  of  penetration  where  the  incident  particle  energy  is  lowest.  The  number  density  of 
displacements  is  found  by  multiplying  the  number  of  primary  knock-ons  by  the  average  number  of 
defects  produced  per  pka,  which  is  calculated  by  averaging  over  the  Rutherford  cross-section-1. 

Figure  1 1  shows  the  damage  profile  in  terms  of  the  number  of  displacements  per  target  atom 
(dpa)  for  23  MeV  alphas  on  tungsten  at  the  dose  used  in  the  preliminary  tests.  As  indicated  above, 
the  damage  peaks  near  the  mean  depth  of  pene  nation.  The  damage  is  about  1  dpa  for  most  of  the 
path,  but  jumps  to  several  hundred  dpa  near  the  end  in  other  words,  each  target  atom  is  knocked 
from  its  lattice  site  several  hundred  times  du’ing  the  activation.  Figure  12  shows  the  distribution 
of  implanted  He  for  the  alpha  activation,  asvunvng  a  gaassian  centered  on  the  mean  depth  of 
penetration  with  the  spread  given  by  Z’^glc;  ■ ' .  tiler  is  bused  on  the  range  straggling  calculations 
of  reference  22.  The  lie  concentration  reaches  a  peak  of  about  5%  or  50, (XX)  ppm  at  the  mean 
depth  of  penetration 

Although  the  detailed  kinetics  have  not  been  examined,  it  is  not  unreasonable  to  conclude  that 
the  problems  encountered  with  the  alpha  activation  were  caused  by  the  combination  of  a  high 
concentration  of  implanted  He  coincident  with  the  highly  damaged  region  at  the  mean  depth  of 
penetration.  The  combination  almost  guarantees  cavity  formation,  which  was  probably  aggravated 
by  the  elevated  temperatures  encountered  during  operation.  The  cavities  may  have  coalesced  into  a 
blister  which  ruptured,  or  the  high  density  of  smaller  discrete  cavities  may  have  decreased  the  heat 
conduction  into  the  cathode  bulk  enough  to  overheat  and  preferentially  melt  the  70  micron  layer 
above  it.  The  jumps  seen  in  the  preliminary  mass  loss  data  could  have  been  caused  by  episodic 
loss  of  gross  amounts  of  material  in  the  molten  layer. 

For  metals,  it  is  recommended  that  die  dose  not  exceed  about  lO16  cm"-  to  avoid  changes  in 
the  target  properties2^.  This  maximum  allowable  dose  places  a  limit  on  the  quantity  of  radioactive 
product  generated  per  unit  area  The  reaction  must  therefore  have  a  high  enough  yield  to  give  a 
sufficient  initial  activity  R  within  this  limit: 

R  =  CAY^t)max  (9) 

where  A  is  the  surface  area  exposed  to  the  beam. 

The  limit  on  allowable  dose  can  be  expressed  as  the  area  which  must  be  irradiated  at  the 
maximum  dose  to  give  an  acceptable  level  of  activity.  In  other  words,  one  can  compensate  for  a 
low  yield  reaction  by  irradiating  a  larger  area.  This  method  of  compensation  obviously  competes 
with  the  need  to  irradiate  small  areas  to  obtain  spatially  resolved  data.  For  the  cathode,  high  spatial 
resolution  is  only  required  in  the  axial  direction  if  axisymmetry  is  assumed.  Therefore,  a 
circumferential  zone  can  be  irradiated  instead  of  a  spot  to  increase  the  area  without  sacrificing  axial 
resolution.  However,  because  of  attenuation  through  the  tungsten  cathode  only  a  certain  fraction 
of  the  total  activity  can  be  observed  perpendicular  to  the  cathode  axis.  For  a  1  cm  diameter  cathode 
this  fraction  is  about  68%. 


vvc-nm'in  m  >  tn  t  a.niyvp.iiiM  ftiTTKl  iFPi 


The  depth  ot  activation  and  the  shape  of  the  radioisotope  density  profile  are  dependent  on  the 


7 


beam  energy  and  die  shape  of  the  cross- section  curve  at  and  below  the  beam  energv.  The 
following  requirements  form  cntena  b\r  mast  lx-  satisfied,  bv  proper  hoi:  of  in>.  idem  p.nw  j. 
energy  and  reaction 


a.  Actuation  Depth 

The  activated  layer  should  be  sufficiently  deep  that  data  can  be  collected  for  a  reasonable- 
period  of  operating  time  and  so  that  the  craters  characteristic  of  arc  damage  do  not  puncture  the 
layer.  These  considerations  must  be  balanced  against  the  required  sensitivity  of  the  technique  Tb 
fractional  change  in  activity  level  is  roughly  the  same  as  the  fraction  ot  the  activated  layer  removed 
Therefore,  if  it  is  possibh  to  detect  a  1%  change  ui  activity,  for  instance,  then  the  sensitivity  is 
limited  to  about  1*T  of  the  activated  depth  and  is  inversely  proportional  the  thickness  of  the  laver 


b.  Activation  Profile 

The  density  profile  of  the  radioactive  (racer  is  important  because  of  its  effect  or  the 
interpretation  of  the  data  The  measured  depth  calibration  curves  assume  mass  loss  in  the  form  1 
thin  uniform  layers  As  explained  above,  a  different  depth  calibration  is  required  d  the  geomtr 
material  loss  differs  from  that  assumed  and  significant  biases  can  he  introduced  by  using  an 
inappropriate  depth  calibration.  Ihe  ratio  o!  the  actual  mass  loss  to  that  cak mated  using,  the  thin 
dice  depth  brat  ion  i- 

Hr)  ;  dV  i  i  ■  A  re  1 

-■  ’  i 

where  T|"  is  the  thickness  of  ,t  -,ln  c  giving  a1'  equivalent  activity  chance 

jv  na  d v  -r  a}  na(r|)  drj  M: 

For  instance,  re  mm  in>:  mng  re-;  u  small  1  .“tuispherical  craters  rather  than  uniform  slices  from 
cathode  av tivated  with  Me-V  <r  pani.de-  produces  the  bias  plotted  in  Figure  i  f  Imfonuw  '■  • 
as  the  plot  show  s,  the  bias  is  a  function  of  starting  depth,  surface  damage  sc  ale,  and  surface 
damage  geometry  For  t're  largest  crater  diameter,  the  underestimation  in  mass  loss  occurs  whet: 
the  crater  punches  through  the  activated  layer  For  the  smaller  craters,  the  effect  is  due  to  the 
non  uniform  distribution  of  the  radioivompe.  As  equation  Ml)  shows,  if  the  radioisotope  density 
were  uniform,  the  volumes  -f  the  actual  ri.m-.igc  H  ‘lie  eqtr'  alen' :  Fee  wauld  tie  equal,  rend<  •  ■ 
the  data  interpretation  independent  of  gem.  '■>. 

The  allowable  deviation  from  a  uniform  density  profile  depends  on  the  desired  accuracy  and 
is  difficult  to  define  precisely,  suite  trie  bias  depends  inherently  on  th”  expected  geometry  of 
material  removal.  Calculations  can  be  performed  to  compare  specific  activation  schemes,  howevci 
In  addition,  non -uniform  erosion  can  lie  audit'd  with  an  activation  scheme  which  produces  two 
radioisotopes  wit!:  d.ft  rent  depth  distrih  t-orre i  otnparison  of  ’he  .  lunges  m  activity  of  th  ■ 
two  sour,  es  can  tv-ea  «>•<’  c-  are  ••  preu  geometry  of  rna-s  loss,  relaxing  to  some  extent  the 
requiremen*  for  prntil*  ui  .-m.iv  and  -coding  an  ever  .rx-r  detailed  picnire  of  the  erosion 
process. 

1- Usable  RiaaipfllMdluilri 

The  radioisotopes  produted  by  the  "  ;u tion  should  meet  these  requirement. 

a.  Half-life. 

The  combined  effect  of  natural  rlccav  and  nuts1-  loss  must  not  reduce  the  activity  below 
dctectibie  limits  before  the  c'r-ettncm  u-  c  pietcd  Tbi-  obviously  d -pends  the  initial  acto  it, 
me  half -ii*e  of  the  i  oto;v,  •> mass  '• *ss  eluting  the  experiment,  these  parameters  can  be 


examined  in  detail  for  a  particular  experiment  and  desired  accuracy,  keeping  in  mind  that  low 
activity  levels  can  be  somewhat  compensated  for  by  longer  counting  times,  but  it  is  a  good  rule  of 
thumb  not  to  plan  an  experiment  lasting  longer  than  three  half  li\ cs. 


b.  Spectral  Separation. 

For  the  least-squares  method  of  spectrum  analyst  to  resolve  spectra  wiJi  multiple  sources, 
the  response  functions  of  those  sources  must  be  linearly  independent.  So,  if  several  isotopes  are 
produced  in  a  single  component,  or  if  several  activated  components  contribute  to  the  spectrum, 
their  individual  spectra  must  be  sufficiently  different  to  allow  precise  resolution. 

c.  Measurable  Response  Functions. 

Another  requirement  dictated  by  the  least  squares  method  is  that  it  be  possible  to  obtain 
single-source  response  functions,  which  is  easily  satisfied  for  activation  schemes  which  yield  onl> 
one  gamma  emitter.  For  those  activations  which  produce  more  than  one  isotope,  however,  there 
must  be  some  method  of  measuring  the  individual  spectra.  For  isotopes  with  different  half-lives 
or  different  depth  profiles  it  may  be  possible  to  i.^olatc  one  isotope  by  allowing  the  shorter  lived 
ones  to  decay  sufficiently  or  by  lapping  off  tnose  with  shallower  profiles.  In  addition,  a  number  of 
pure  sources  are  available  commercially. 

THE  IMPROVED  ACTIVATION  PROCEDURE 

Table  3  compares  the  tungsten  activation  scheme  used  in  the  preliminary  tests  with  the 
deuteron  activation  chosen  for  subsequent  experiments  on  the  basis  of  the  criteria  given  above. 

The  results  of  the  preliminary  tests  with  the  a-activation  demonstrate  the  consequences  of 
exceeding  the  dose  limit  for  avoiding  structural  damage.  The  new  deuteron  activation  produces 
several  isotopes  with  a  much  larger  yield,  affording  extremely  good  spatial  resolution  within  the 
allowable  dose  limits.  Figure  14  shows  that  the  deuteron  activation  produces  less  than  10' 3  dpa 
for  most  of  the  depth  and  peaks  at  only  about  1  dpa  below  the  activated  layer,  permitting  the 
conclusion  that  radiation  damage  and  effects  will  not  affect  the  results  of  the  measurement 

The  depth  of  activation  for  the  184Re  is  about  100  microns,  much  thicker  than  that  of  the  it 
reactions.  This  allows  the  benefits  of  thicker  layers,  but  sacrifices  some  depth  resolution.  To 
regain  that  resolution  however,  the  activated  layer  can  be  compressed  by  irradiating  at  an  angle  ' 
the  surface.  Figure  15  displays  the  depth  calibration  curve  for  this  reaction  measured  on  a  targei 
irradiated  at  an  angle  of  75°  from  the  normal,  compressing  the  layer  by  a  factor  of  4. 

The  linearity  of  the  depth  calibration  curve  reflects  the  more  uniform  density  of  radioactive 
tracer  associated  with  this  activation.  The  ratio  of  actual  mass  loss  to  measured  mass  loss  for  this 
profile  is  plotted  in  Figure  16,  demonstrating  much  less  non-uniformity  bias  than  the  a  activation. 

It  was  hoped  that  the  *^2Ta  produced  in  this  activation  could  be  used  as  a  second  layer  to  study  the 
surface  damage  geometry,  but  the  reaction  yield  is  too  low  to  provide  usable  levels  of  activity.  It 
may  still  be  possible  to  use  the  182Re  since  it  has  a  sufficiently  high  yield,  but  its  short  half-life 
severely  limits  its  useful  lifetime. 

The  deuteron  activation  produces  a  number  of  other  intense  isotopes,  but  most  of  these  are 
very  short-lived  and  decay  to  negligible  levels  within  days.  About  3  weeks  after  activation  the 
target  can  be  assumed  to  be  a  pure384Re  source,  which  allows  simple  generation  of  a  reference 
response  function.  The  i82Re  spectrum  is  sufficiently  different  from  die  184Re  spectrum  to  allow 
good  resolution  with  the  linear  least  squares  approach.  The  !82Re  response  function  can  be 
generated  by  subtracting  the  184Re  response  function  from  the  composite  spectrum. 

The  polynomial  fit  used  to  interpolate  the  eroded  depth  from  the  depth  calibration  data  yields 
values  with  a  standard  error  of  0  2  micrometers,  which  represents  an  approximation  of  the 
statistical  error  inherent  in  this  technique.  In  addition,  a  systematic  error  is  introduced  by  using  the 
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thin- slice  depth  calibration  to  inter;-.'?!  tV  data  1  lowevei .  fie  ire  I  b  demonstrates  that  th;%  en  '  • 
less  than  a  feu  percent  for  most  ot  the  dept!-  range  for  the  smaller  craters  and  is  a  sipr.ific.in!  ’  •- 
only  for  those  craters  which  punch  thnnigh  'he  layer  In  add".ion.  It  must  be  recognized  the  e 
values  can  he  reduced  even  h’rthcr  if  an  activated  depth  of  more  than  ">5  micrometers  is  used 
Finally,  the  vise  of  the  ,x*Re  tn  conjunction  with  the  'wRe  mas  eliminate  the  uncertain  tv 
introduced  by  the  nonuniform  geometry  of  cat*vxie  surface  .iamage 

R  hUF  APPLICATION  OK  THE  REFINED  SLA  1  FCHNIQI'F 

Although  it  is  important  to  understand  the  en>sion  mechanisms  dominating  in  the  microspot 
emission  mode,  the  unacceptably  high  eros;on  rate  associated  with  this  emission  mechanism 
mandates  exploration  of  the  potentially  more  benign  thermionically  emitting  cathixJe  Two 
independent  paths  will  be  followed  to  investigate  this  regime  of  cathode  operation --«miuh!ioe 
the  hot  cathode  of  a  steady  state  device  using  a  multi-megawatt  pulsed  MPI)  thruster  and  actual 
time-resolved  erosion  measurements  on  a  low-power  steady  state  thruster 

PULSED  MPD  THRUSTER  EROSION  MEASUREMENTS 

A  half-scale  benchmark  thruster2  **  similar  to  that  shown  in  Figure  17  will  lx-  used  for  thi 
series  of  tests.  After  obtaining  a  more  reliable  set  of  erosion  data  from  this  pulsed  device  with 
bulk  cathode  surface  at  nlxc  roonvtempcriture  to  characterize  the  m:cn»sp>t  emission  mode.  n 
sequence  of  testing  will  performed  with  the  cathode  externally  heated  prior  to  the  discharge  M 
an  inductive  coil  This  precise  control  over  the  bulk  surface  temperature  will  allow  examining  pic 
transition  from  a  cold  cathode  dominated  by  microspot  emission  to  an  incandescent  •  athod?  on 
which  current  continuity  is  satis Hed  by  th/  muonic  emission 

TIME-RESOLVED  STEADY  STATE  THRUSTER  EROSION  TESTS 

lire  real-time,  time  resolved  erosion  measurement  technique  will  be  used  to  study  the  inm 
high-erosiori  start  up  phase  and  the  less  destructive  steady-state  operation  on  the  coaxial  dev  it ; 
pictured  in  Figure  18,  which  is  an  MPD  type-  thruster  operated  at  power  levels  of  up  to  30  kW  .c 
current  levels  from  500  to  1000  A  for  periods  of  2  minutes  and  less26.  The  cathode  in  this  thnotc 
experiences  current  densities  and  surface  temperatures  similar  tn  those  expected  in  high  power 
steady  state  MPI)  thrusters  Measurement  t-t  the  local  cathode  tempemture  in 
conjunction  with  the  mass  loss  measure*"  cro  will  allow  cotupanvin  with  erosion  data  taken  i  — 
the  external!  v  heated  cs’Ikoc  of  ’he  p  .,;!••.«  >  MPD  thruMet 

CONOAfSlONS 

Preliminary  tests  of  the  SLA  technique  confirm  that  pulsed  thrusters  for  which  bulk  cathode 
temperatures  remain  below  that  required  for  significant  thermionic  emission  are  subject  to  severe 
erosion  associated  with  the  miens  me  rn  su  m:  mechanism,  which  renders  them  useless  as 
practical  tht usu  ry.  Pub.  d  dev  iv.  toe:  il  foi  simulating  tbc  opener  -u  of  steady  state  thru- 
however,  and  remain  an  integral  p.v?  of  :‘x  experimental  program  Anomalous  jumps  in  the 
erosion  data  and  the  blistered  appeiinmc"  •>  the  activated  spot  can  be  at'iitntted  to  structural 
radiation  damage  caused  by  the  high  doses  of  a- particles  used  in  the  activation 

The  criteria  outlined  above  are  desmned  to  addtess  the  problem  of  radiation  damage  and 
several  deficiencies  in  the  data  red  at  riot:  T,;  Jcu'.eroa  activation  chosen  on  the  basis  of  these 
criteria  appears  capable  of  providing  extrertv.  iy  accurate,  non  intrusive  measurements  of  cathode 
erosion.  The  dcuteion  activation  paj.itia-vtv  and  depth  calibration  c  urve  given  above  can  be  u  1 
with  tungsten  surhn.es  activated  in  a  umvcrci'v  or  government  cyclotn*n  facility  or  from  the 
commercial  activation  service  offered  b\  tl  V  K  Harwell  Laboratory  in  Harwell,  England. 
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Lrosion  measurements  can  then  be  performed  with  standard  radiation  measurement  equipment 
similar  to  that  described  in  the  first  section  of  this  paper  and  analyzed  on  a  microcomputer  with  the 
gamma  spectrum  analysis  programs  noted  in  the  references.  The  criteria  outlined  in  the  second 
section  of  the  paper  can  be  used  to  guide  the  development  of  an  activation  scheme  appropriate  for 
other  materials  of  interest.  Activation,  data  collection,  and  interpretation  have  now  reached  a 
sufficient  level  of  sophistication  that  surface  layer  activation  can  become  a  routine  diagnostic 
technique  for  studying  surface  degradation  on  any  material  that  can  be  suitably  activated. 

The  planned  erosion  measurements  on  the  pulsed  thruster  with  the  externally  heated  cathode 
and  the  steady  state  thruster  should  provide  a  sound  data  base  for  directing  and  testing  erosion 
modelling  efforts.  In  addition,  the  time-resolved  erosion  measurement  will  demonstrate  a  new 
application  of  surface  layer  activation  that  can  provide  an  extremely  powerful  method  of  lifetime 
testing  steady  state  thrusters. 
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Table  3: 

Comparison  of  Tungsten  Activations 

Beam 

Beam 

Energy 

(MeV)  Reactions 

Yield 

(ptCi/C) 

Activated  Axial 

Area*  Resolution** 
(cm2)  (Cm) 

Activated 
Depth  t 
(pun) 

Half-life 

(days) 

Major 

y  1  .ines 

(kcV) 

a 

23  182\V(a.n)1850s 
i83w(a,2n)185os 
184w(a3n)185os 
186w(a,5n)1850s 

50 

6.25 

2.93 

26 

93.6 

646 

d 

15  182\v(d,Y)lg4Re 
183w(dji)184Re 
184\V(d,2n)184Re 
186w(d,4n)>84Re 

7000 

0,05 

0.02 

100 

38.0 

-no 

1  J  4- 

895 

903 

182w(d,2p)l82Ta 

184w(d,a)182Ta 

40 

8  75 

3.50 

70 

115.0 

1121,  1189 
1221, 1231 

I82w(d»I82RC 

183w(d,3n)182Re 

184wfd,4n)182Re 

15,000 

0  02 

O.Oi 

L 

2.67 

toy,  229 
256,  351 
1076, 1121 

*  Area  for  1  ptCi  activity  at  a  dose  of  10*6  cm*2 

**  Axial  length  of  circumferential  activated  strip,  assuming  area  from  column  6  and  an  attenuation  of  32^ 
t  For  beam  normal  to  the  surface 
+  Not  yet  measured 
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Figure  2:  Composite  gamma-ray  spectrum. 
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calibration  for  tungsten  activated  with  23 
:les  normal  to  the  surface. 


;  10:  Comparison  of  calculated  (— )  and  measured  ( 
calibration  curves  for  *l-.51  MeV  protons  normal  to 


Figure  12:  Profile 
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Figure  16:  Non-uniform  erosion  bias  for  tungsten  deuteron 
activation. 
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Figure 


